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Abstract. Mathematical justifications are given for several integral and series 
representations of the Dirac delta function which appear in the physics liter- 
ature. These include integrals of products of Airy functions, and of Coulomb 
wave functions; they also include series of products of Laguerre polynomials 
and of spherical harmonics. The methods used are essentially based on the 
asymptotic behavior of these special functions. 

1. Introduction 

The Dirac delta function S{x) has been used in physics well before the theory 
of distributions (generalized functions) was introduced by mathematicians. The 
manner in which physicists used this function was to define it by the equations 

(1.1) 6{x — a) = 0, X a, 

and 

/oo 
(j){x)6{x - a)dx ^ (j>{a), a e E, 

-co 

for any continuous function </'(x) on M. However, mathematically, these two equa- 
tions are inconsistent in the classical sense of a function and an integral, since the 
value of the integral of a function which is zero everywhere except for a finite number 
of points should be zero. There are now two mathematically meaningful approaches 
to help us interpret the delta function given in (1.1) — (1-2). One approach is to 
consider 6a '■= S{x — a) as a continuous linear functional acting on a space of smooth 
functions with rapid decay at ±00, and the action of Sa on a particular function 
(j){x) is given the value (j){a); see [14, p. 141] and [ , p. 77]. The other approach is 
to find a sequence of functions S„{x — a) such that 

/oo 
Sn{x — a)(f){x)dx = 0(a), a G M; 

-00 

see [7, p. 55] and [i. , p. 17]. Such a sequence is called a delta sequence and we write, 
symbolically, 

(1.4) \\m 5n{x - a) ^ 5{x - a), x e M. 
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It seems that the second approach is more acceptable to physicists and apphed 
mathematicians . 

Recently, in the process of preparing some material for the major project "NIST 
Handbook of Mathematical Functions [I-:)]", we encountered some very interesting 
integral and series representations of the delta function which need mathematical 
justification. For instance, in [ , p. 696] the formula 

/•CO 

(1.5) / xtJi,{xt)J^{at)dt — 5{x ~ a), Re > — 1, x > 0, a > 0, 



appears, where Jy{x) is the Bessel function of the first kind, and in [" , Eq.(122)] 
one finds the integral representation 

/>C30 

(1.6) / s{x,l;t)s{a,l;t)dt — S{x — a), a > 0, x > 0, 



where s{x,l;t) is the Coulomb wave function. A recent reference [I V, p. 57] on the 
Airy function Ai(a;) also gives the formula 

/oo 
Ai{t - x)Ai{t - a)dt = 5{x - a). 
'OO 

While physicists may find these representations convenient to use in applications, 
mathematicians would, in general, feel uneasy or even disturbed to see these for- 
mulas being used since the integrals in (1.5) — (1.7) are all divergent. Thus, it 
would seem meaningful and necessary to give a mathematical justification for these 
representations, and this is exactly the purpose of the present paper. 

There are also some series representations for the delta function. These include 
the following: 



(1.8) y [k+-\Pk{x)Pk{a)^5{x-a), 



— {x+a)/2 °° 1 

(1-9) ^^i7,(a;)i7,(a) = 5(x-a), 

fe=0 



and 

oo 

(1.10) e-("+'')/2 ^ Lk{x)Lk{a) = S{x - a), 

fe=0 

where P}~{x), Hk{x) and Lk{x) are, respectively, the Legendre, Hermite and La- 
guerre polynomials. Equations (1.8) — (1.10) are special cases of an equation in 
Morse and Feshbach [I J , p. 729]. Another series representation is given in [-S, p. 792]; 
that is, 

oo k -J 

(1-11) t^mt^k '"^^1 

= J(cOs6'i — COS02)<5((/>1 — 4>2), 

where the functions Yki{d,4>) are the spherical harmonics (see [3, p. 788]) and the 
asterisk "*" denotes complex conjugate. 

The orthogonal polynomials in (1.8) — (1.10) and the orthogonal function in (1.11) 
are the eigenfunctions corresponding to the eigenvalues (discrete spectrum) of some 
differential operators. Likewise, the special functions in (1.5) — (1.7) can be regarded 
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as the eigenfunctions associated with the continuous spectrum of corresponding 
differential operators. The proofs of the representations in (1.8) — (1.11) turn out 
to be much simpler than the proofs of those in (1.5) — (1-7). Indeed, we shall show 
that the results in (1.8) — (1.11) all follow from expansion theorems in orthogonal 
polynomials, whereas for the representations in (1.5) — (1.7) we need to provide 
some new arguments. 

2. A GENERALIZED RIEMANN-LEBESGUE LEMMA 

There are already several delta sequences in the literature. For instance, we have 

(2.1) 5„{x - a) = ^e-"(—^\ 

Ti 1 

(2.2) 5^{x-a)^--—^ ^, 

and 

^ , ,1 sin n(x — a) 

2.3 5r,{x-a)^ ^ '-■ 

TT X — a 

see [5, pp. 35-38] and [8, pp. 5-13]. To verify whether a given sequence of functions 
is a delta sequence, one can apply the criteria given in [5, p. 34]. If the function (^{x) 
in (1.3) is only piecewise continuous in K, then this equation becomes 

/■°° 1 

(2.4) lim / 5Jx-a)<j)(x)dx= -[(l){a+) + <j){a^% a e M; 

n^ooj_^ 2 

see [8, p. 16]. 

For convenience in our later argument, we also state and prove the following 
result. 

Lemma. (A generalized Riemann-Lebesgue lemma). Let g{x,R) be a continuous 
function of x £ (A, B) and uniformly bounded for R > 0. If 

(2.5) lim / g{x, R)dx = 

R-^°^ J A' 

for any A' and B' with A < A' < B' < B , then 

(2.6) lim / 'i{;{x)g{x,R)dx = 

for any integrable function ^p{x) on the finite interval {A, B). If A = and B = +oo, 
or if A = —oo and B = +oo, then (2.6) holds for any absolutely integrable function 
ip{x) on the infinite interval {A,B). 

Proof. First, from (2.5) it is easy to see that (2.6) holds for step functions. Now, 
let 'ip{x) be an integrable function on {A,B). For any e > 0, we can always find a 
step function s{x) such that 

(2.7) £\^x)-s{x)\dx<^, 

where K = max{|5(a;, R)\ : x G M and R > 0}. Choose i?o > so that 



(2.8) 



s{x)g{x, R)dx 



< I for all R > Rq. 



4 



Y. T. LI AND R. WONG 



Write 



I '^{x)g{x,R)dx = I ['ij:{x) — s{x)\g{x,R)d: 

J A J A 

On account of (2.7) and (2.8), we have 

4>{x)g{x, R)dx 



H{x)g{x, R)dx. 



< e 

for all R > i?o- Since e is arbitrary, this proves (2.6) when A and B are finite. 

If the interval of integration is infinite, and if ip{x) is absolutely integrablc there, 
then we can choose finite numbers A and B such that the integral / ip{x)g{x, R)dx 
outside the interval (A, B) is small since g{x, R) is uniformly bounded. On the finite 
interval {A,B), we can apply the result just established. □ 

3. Bessel function 
The Bessel function Ju{xt) is a solution of the differential equation 



(3.1) 



d^fdy 



dt \ dt 



+ [xH--]y = 0. 



With X replaced by a, one obtains a corresponding equation for J^{at). Multiplying 
equation (3.1) by Jv{at) and the corresponding equation for Jv{at) by Jv{xt), and 
subtracting the two resulting equations, leads to 



d 



(3.2) 



[x^ — a^)tJ„{at)J„{xt) = —{atJi,{xt)Jl{at) 



xtJv{at)Jl{xt)}. 



Put 

(3.3) 5R{x,a)=x 
An integration of (3.2) gives 



J„{xt)J^{at)tdt. 



(3.4) 



5R{x,a) = 



at Ju{xt)Jl{at) — xt Jv{at)Jl{xt) 



Prom the ascending power series representation 

(-1)" 



(3.5) 



^ ' 71=0 



T{n + u + l)n\ \2 



2n 



it can be shown that the leading terms in the series expansions of at Jv{xt)Jl{at) 
and xt J^{at)Jl{xt) cancel out. Thus, the right-hand side of (3.4) vanishes at the 
lower limit when v > —1, and we obtain 



(3.6) 



6R{x,a) 



;[aRJ^{xR)Jl{aR) - xRJ^{aR)Jl{xR)]. 



Theorem 1. For a > 0, i' > —1 and any piecewise continuously differentiable 

function 4>(x) on (0,00), we have 

(3.7) lim / (f)(x)x( [ J„(xt)J^(at)tdt^dx=h(l)(a+) + (j)(a~)], 
provided that 
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(i) x~^(j){x)\dx converges; 

(ii) x^(j){x)\dx converges when v > —\, or 
(ii') x^^^\(f){x)\dx converges when —1 < v 

Proof. In view of the asymptotic formulas 



(3.8) 

and 

(3.9) 



Jv{xt) = 



Jlixt) 



2 

TTXt 



2 
TTXt 



cos( - Y~ 4 



sin I xt 



+ si{x,t) 



■ £2ix,t) 



where ej{x,t) = 0(l/t) as t ^ oo uniformly for x > 5 > and j = 1,2, there are 
constants Mi > and M2 > such that 



\J^{aR)\ < Mii?-5, 



|4(ai?)| < MiR-i 



(3.10) 
and 

(3.11) \JuixR)\ < M2X-^R-^, \JUxR)\ < M2X-^R-i 

for a; > 1 and R> 1. From (3.6), it follows that 



\6R{x,a)\<MiM2-. 



.3/2 



+ 



1 



;V2 ; ■ 



Hence, for b > max{a, 1}, we have 

\6R{x,a)(f){x)\dx < M3 I x~i\(l){x)\dx, 



where M3 = MiM2{a + l)&^/(6^ — a^). Since the last integral is convergent by 
condition (i), for any e > there exists a number c>b such that 



(3.12) 



j \5R{x,a)(j){x)\dx < 



Let < p < min{a, 1}. To estimate the integral of 5fi{x,a)4>{x) on the interval 
(0,/o), we divide our discussion into two cases: (i) v> — i, and (ii) —l<v<—^. 
In the first case, we have from (3.5) a positive constant M4 such that 

(3.13a) \Ju{xR)\ < M'^{xRY < Mi{xR)-i 

and 

(3.13b) \Jl{xR)\ < MiixR)"-^ < M'^x-iR-i < M'^x'^R'^^ 

for < xR < 1 and R>1. Prom (3.8) and (3.9), we also have 

(3.14) \Ju{xR)\ < M'lx-^R-^, \Jl{xR)\ < M'^x-iR-i 

for xR > 1. Coupling (3.13) and (3.14), we obtain 

\MxR)\ < M4X^iR~i, \Jl{xR)\ < MiX^iR-^ 

for < a; < 1 and i? > 1, with M4 = max{M^, M'^}. Thus, 

1 

X2 



|<5fl(x,a)| <MiM4(a + l) 



a? — p^^ 



< a; < 1, 
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and 

(3.15) r \6R{x,a)<l){x)\dx < M5 f xi\(j){x)\dx, 

Jo Jo 

where M5 = MiM4(a + l)/{a'^ - p^). 

In the second case, there are constants Mg > and Mg > such that 

\Ju{xR)\ < M^ixR)" < M^x^R-K 

\JUxR)\ < M^ixR)"-'^ < M^x^-^R-^ 

for < xR < 1 and -R > 1, and 

\MxR)\ < M^x-^R-^ < M^x^R-^ 

\JUxR)\ < Ml^x-iR-i < M^V-ii?-5 

for a;E > 1 and ii > 1. With Me = max{M^, M^}, it follows that 

\MxR)\ < Mex^R-^ and \JUxR)\ < Mqx^-^R-^^. 

Therefore, 



1 

a" — 



\5R{x,a)\ < MiM(,^—^{a + l)x"+^ 



and 

fp fP 

(3.16) / \SR{x,a)(f){x)\dx < Mr x''+^\(t){x)\dx, 
Jo Jo 

where My = MiM(i{a. + l)/{a'^ ^ P^)- On account of conditions (ii) and (ii'), for any 
£ > there is a constant < rf < p such that 

(3.17) \5R{x,a)ct>{x)\dx < -. 

For d < X < c, & combination of (3.6), (3.8) and (3.9) gives 

2 X { fx 
dR{x,a) = 2 2 1 \/ ~ 



(3.18) 



^ cos ((x, R) sin ({a, R) + e(x, a;R) > , 
X J 



where c and d arc given in (3.12) and (3.17), respectively, <^{x,R) = xR— ^ — f 
and e{x,a;R) = 0{1/R) as i? — > 00 uniformly for x € (rf, c). Note that e{x,a;R) 
is continuously differentiable, e{x,a; R)/{x — a) is continuous in x and uniformly 
bounded in x and R, and (f){x) is piecewise continuous in {d,c). Hence, 

(3.19) lim r ^K c^(x)dx = 0. 

J^ x^ - 0? 

By the Riemann-Lebesgue lemma, we also have 
hm /Y/-l)^, 



(3.20) lim r (J^--^''^^^^^2^x4,{x)dx = ^ 



and 
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A combination of the results in (3.18) — (3.21) yields 

lim / dii(x, a)(p[x)ax — lim — / ^ ^ x<p(x)ax 



R—>OD J ^ i?.— ^oo TT 

f^svuRix — a) 2x 

= lim / r — • (p[x)ax 

if^ooj^ 7r[x ~ a) x + a 

On account of Jordan's theorem on the Dirichlet kernel p. 473], we conclude 

f 1 
(3.22) lim / dR{x,a)(t){x)dx^ -[(l){a') + (f>ia+)]. 

Since the number e in (3.12) and (3.17) is arbitrary, it follows from (3.22) that 

f°° 1 _ 

lim / dR{x,a)(t){x)dx = -[(j){a ) + 0(a+)], 

which is equivalent to (3.7). □ 

4. Coulomb wave function 
The Coulomb wave function s(a;, /; r) is a solution of the Coulomb wave equation 

(4-1) d^+^+u- 

that satisfies the initial conditions 

(4.2) s{x,l;0) = s'ixj;0) ^0 
and has the asymptotic behavior 

(4.3) s{x, I; r) = — ^a;^i [sin^(a::, /; r) + £i(x, I; r)], 

(4.4) s'{x,l;r) = -^xi[cosC{x, l;r) +e2{x,l;r)], 

where Sj{x, I; r) = 0{l/r) as r — > oo uniformly for x > d > 0, j — 1,2, and 

(4.5) C{x,l;r)^ kr +\\n {2kr) - — + arg T | ^ + 1 - | 

k 2 \ k 

with k — s/x; see [l(i, p. 236]. For a > and a; > 0, we define 

(4.6) 3ji(x,d)= / s{x,l;r)s{a,l;r)dr. 

Jo 

From (4.1) and (4.2), one can show as in Sec. 3 that 

s{x,l;R)s'{a,l;R) - s'{x,l;R)s{a,l;R) 

(4.7) dR(x,a)^ . 

X ^ a 

Theorem 2. For any a > and any piecewise continuously differentiable function 
4i{x) on (0, cx)), we have 

(4.8) lim / (j){x) I s{x,l;r)s{a,l;r)drdx=-[(l){a^) + (j){a~)], 
R-^°oJo Jo 2 

provided that the integrals 



(4.9) / X ^\4'{x)\dx and / x i\4){x)\dx 

are convergent. 
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Proof. From formulas (4.3) and (4.4), for any b > max{a, 1} there exists a number 
Ml > such that 

\s{x,l;R)s'{a,l;R)\ < Mix^i, 
\s'{x,l;R)s{a,l;R)\ < Mix^ 
for i? > 1 and x > b. Hence, it follows from (4.7) that 

\SR{x,a)c^{x)\dx<2Mi H ^^^^^dx < M2 x"! |0(x)|dx, 
Jb x-a 7b 

where M2 — 2bMi/{b — a). By hypothesis, the last integral is convergent; so for 
any e > there is a number c > b such that 

(4.10) / \dR{x,a)(t>ix)\dx < - for alii? >1. 
To prove that there exists a number d > such that 

(4.11) / \bB.{x,a)(l>{x)\dx < - foraUi?>l, 
Jo 2 

we first need to demonstrate that 

(4.12) fc5s(fc^;;r) = 0(1) as r ^- cx), 

(4.13) s'(fc2,?;r) = 0(1) as r ^ cxd, 

uniformly for all sufRciently small fc > 0. (Recall: fc = \fx.) This can be done 
by considering two separate cases: (i) fcr — >■ cx), and (ii) kr bounded. In case (i), 
we first make the change of variable p = kr and set a;(p) = y{p/k) = y{r) so that 
equation (4.1) becomes 

+ fi + A _ fill 1) 

dp^ \ kp pi- 



and then apply the Liouville- Green transformation given in [12, p. 196] with /(p) = 
1 + (2/fcp) and g{p) = 1(1 + l)/p^. The result is that equation (4.14) has a solution 
uJi{k; p) such that 

(4.15) wi(fc; p) ^1 + sin| /" (^1 + dp\ as p 00 

and 



kp J VJ \ 



(4.16) |.i(fc; P) - (1 + cos{ / (1 + ^) ^ dp] as p ^ 00. 

With a suitable choice of the integration constant and for fixed fc > 0, we have 

1 + ^^ dp = C*(fc,^;/o) + oQ^ asp^-oo, 
where 

1 , Ztt /, i 



C(fc,Z;p) =p+^ln2p - - +argr( / + 1- - 
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which is exactly equal to the function C,{k'^,l\r) given in (4.5). For fixed A: > 0, we 
can compare the behavior of s{k'^,l;r) given in (4.3) and that given in (4.15). The 
conclusion is 



(4.17) sik',l;r) 
from which we also obtain 



1 

/ Trfc 



(4.18) 



wi(fc;p), 



k duJi 
TT dp 



Since a;i(fc; p) = 0(1) and id'i{k: p) = 0{k^i) for all small k and large p on account 
of (4.15) and (4.16), the order estimates in (4.12) and (4.13) are established. 
In case (ii), we first recall the function 



{^|k) 



r(2Z + 2) 



M,/,^,+ i(-2zfcr), 



where Mi^^x{z) is a Whittaker function; see Seaton [10, eqs. (14) & (22)]. This 
function is related to the Coulomb wave function s(fc^, I; r) via 



s{k\l-r) 



2(1 - e-2Vfe) 



where A{k'^,l) is a polynomial of degree / in /c^; see [Ki, eq.(114)]. In view of the 
convergent expansion [ , § 7, eq.(16)] 



where the pl^^'' (z) are polynomials in , we have for bounded kr 



.,.(.) = r(2A + i)2^V+^^p(^^)(z) ^;^^^^) , 



n=0 



f{k^,l;r) = C{kr)ri cos( 2^2^- |(2/ + 1) - ^tt 



as r — cx), where C{kr) is a polynomial of kr. Hence, 

A{k'^,l) 



O 



1 

HV4 



(4.19) 



and 



(4.20) 



s{k\l\r) = C*ikr)k~^ 



2(1 - e-2Vfe) 



cos( 2V2r - -{21 + 1) 



O 



1 

HV4 



s'{k\l]r) = -C{kr)r- 



(1 - e-2Vfc) 



sin( 2V2r- -(2/ + 1) 



1 



where C*{kr) = {kr)^C{kr), again proving (4.12) and (4.13). 

From (4.12) and (4.13), it follows that there are constants < p < min{a, 1}, 
i?o > and iVi > such that 



\xis{x,l\r)\ < Ni 



and 



\s'{xJ;r)\<N^ 
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for all r > Rq and < a; < p. Furthermore, by (4.7), 

\SR{x,a)\ < ^^{x-i+a-i) < 2N!j^ 
\x — a\ \x — a\ 

for < X < p and R > Rq, and 

fP fP _^ 

/ \Si{{x,a)4>{x)\dx < N2 / X ■i\4>{x)\dx, 
Jo Jo 

where N2 = 2Nf/{a — p). By hypothesis, the last integral is convergent, thus 
establishing (4.11). 

Let us now consider the case when x lies in the interval {d, c). Prom (4.3), (4.4) 
and (4.7), we have 

tt{x — a)Sii{x, a) — i — ] sin ({x, l; R) cos ({a, I; R) 



(4.21) 

COS ({x, I; R) sin ({a, I; R) + e{x, a; R), 



where e{x, a; R)/{x — a) is continuous in (0, 00) and e{x, a; R) = 0{l/R) as R ^ 00 
uniformly for x G {d,c). As a consequence, we obtain 

(4.22) lim r (j){x)dx = 0. 

X — a 

For any A and B satisfying d < A < B < c, a,n integration by parts yields 



cosC(A;2,/;i?) 



A 



where Cfe(fc^, I', R) denotes the derivative of C(^^) I', R) with respect to k. Prom (4.5), 
it is readily seen that for d < k'^ < c, 

a{kM;R)^ 00 and |(^^(^)-^0 

as — >■ 00. Hence, 

(4.23) lim / smC{k^,l;R)dk = 0. 
Similarly, we also have 

(4.24) lim / cosC(fc^Z;i^)dfc = 0. 

Let ?7 > be an arbitrary number such that d<a — ri<a + ri<c. A combination 
of (4.23), (4.24) and the generalized Riemann-Lebesgue lemma gives 

(4.25) hm ( P V ](^\ ^^smC{x,l;R)dx = 
and 



\ .Id Ja+riJ \ ^ J X - a 



(4.26) lim ( / +/ W-Y ^^cosC{x,l;R)dx = 0. 



/ + / 

J d J a 
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By the same reasoning, we have 

fa+T] 



(4.27) 

and 

(4.28) 



Urn I 



"^^^^ smC,{x,l;R)dx = Q 



X — a 



hm 











' a—Tj 





cPix) 



X — a 



cos ({x, I; R)dx = 0. 



Prom (4.21), (4.22) and (4.25) - (4.28), it follows that 
lim / 5ji{x,a)(j){x)dx 

ra+ri 

= lim / 

' J a— 7} 



(4.29) 



sm{C{x,l;R)-C{a,l;R)} 



= lim 

R—^oo 



X — a 



(j){x)dx 



j-a ra+ri 
J a—T) Ja 



sin{C(x, i?) - C,{a,l;R)} 



^{x)dx. 



Since (j){x) is piecewise continuously difFerentiable in (0,00), it is continuously 
differentiate in (a — 77, a) for sufficiently small rj > 0, and (j){x) — (j){a~)/{x — a) is 
integrable on (a — r],a). By (4.23), (4.24) and the generalized Riemann-Lebesgue 
lemma, we have 



(4.30) 



lim 

R—^oo 



(x) - 0(a ) 



a— 7] 



X — a 



sin{C(a;, I; R) - C{a, I; R)}dx = 0, 



or equivalently 

sm{C{xJ;R) -CiaJ]R)} 



lim 



^ a—T) 



7r(a; — o) 



(f){x)dx 



r 

= 4>{a ) lim / 



" sm{C{x,l;R)-(i{a,l;R)} 



tt{x — a) 



dx. 



To obtain the value of the limit on the right-hand side of the last equation, we 
shall use the Cauchy residue theorem. Let 



(4.31) 

so that 
(4.32) 



C(x, R) = y/xR 
1 



1 



\nR, 



l-K 



e{x,l) = -= ln(2\/S) - — + argr / + 1 - 



C{xJ;R) = C{x;R)+eix,l). 



Furthermore, let T denote the positively oriented closed contour depicted in Figure 
1 below. It consists of a horizontal line segment T3, two vertical line segments r2 
and a quarter-circle S centered at z = a with radius r, and the interval Fi on 
the positive real-axis. The entire region bounded by F lies in the first quadrant 
e C : Re ^; > and Im 2; > 0}. Consider the complex- value function 



(4.33) 



Fr{z) = 



gi[C(2,iJ)-f(a,fl)] 



z — a 
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Since C(-z, R) is analytic in the right half-plane, by Cauchy's theorem 



(4.34) 



FR{z)dz = 0. 



By Cauchy's residue theorem, we also have 
(4.35) 

where r is the radius of the quarter-circle S 



lim / Fii(z)dz = —in 11^ 




a-y\ 



Figure 1. Contour F. 



For z S F3, we write z = u + iD. Since ({a, R) is real, a simple estimation gives 



(4.36) 



\Fr{z)\ < le-i-C(^,fl)^ 



From (4.31), it is readily seen that limfl_>oo Im C(z,i?) — +00 uniformly for z G F3. 
Hence, lim^_j.oo e~^'° = uniformly for z e F3. From (4.36), it follows that 



(4.37) 



lim 

-R— >-oo 



FR{z)dz = 0. 



For 2; e F4, we write z ~ a — rj + iv. Clearly 



(4.38) 



\Fr{z)\ < ie-i-C(^.«), 



z e r4. 



Let CTi > be any small number. From (4.31), we have lim/j_>.oo Im C(z,i?) = +00 
uniformly for z G r4 and Im z — v > <7i. Thus, lim/j_>oo e~'™ Ciz,R) — q uniformly 
for z G r4 n {z : Im z > ai} and 



lim 

i?.— >oo Jy 



g -imC(.,fl)^^^ lim / e-i'^^^^'^^dz. 



a— 7] 



From (4.31), it also follows that there is a constant M5 > such that e~^™ C(z,-R) < 
M5. Hence 

lim / e-^'^^^^'-'^Uz <Ahai. 
Since cti can be arbitrarily small, we obtain 



(4.39) 



lim 

-R— foo 
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Coupling (4.38) and (4.39) gives 

(4.40) lim / Fn{z)dz = 0. 
In a similar manner, one can establish 

(4.41) lim / FR{z)dz = 0. 

By a combination of (4.34), (4.35), (4.37) and (4.40) - (4.41), we obtain 



lim r ^MC(.>^)-CK^)},, 



(4.42) ""^^^Ja-r, 7r{x-a) 



-Im I lim lim / FR(z)dz\ = -. 

TT [fl-)-oor-i.O+ J 2 



Using (4.42), wc arc now ready to handle the limit on the right-hand side of the 
equation following (4.30). 

For simplicity, let us write 6{x) = 9{x, I); cf. (4.31). By the mean- value theorem, 
cos[^(x) — 9{a)] = l + 0{{x — a)^} for x near a. Hence, by the generalized Riemann- 
Lebesgue lemma, we have from (4.42) 

r MCjx, R) - ga, R)} cos[g(rr) - e{a)] ^ 1 
R^ooJ^_^ n{x-a) 2' 

and also from (4.23) and (4.24) 

(4.44) lim /" cos{ax,R)--aa,R)}sm[eix)-e{a)] ^^_^_ 
R^oo n{x - a) 

Upon using an addition formula, it follows from (4.43) and (4.44) that 

(4.45) lim r Bin{C(x,/;i^)-C(a,/;i^)},^,l_ 
R^oo 7r(a; — o) 2 

Couphng (4.30) and (4.45), we obtain 

/. T r sm{C(x,l;R)-ga,l;R)} ^ , , 

(4.46) lim / — ' ' ' — V <b(x)dx = -Ma'). 

In a similar manner, we also have 

,• r^''sin{C(:r,l;i?)-C(a,«;i?)}^, ,j 1^,+, 
4.47 hm / ' I ' Y " ^{x)dx = -<t>{a+). 

A combination of (4.29), (4.46) and (4.47) yields 

1 

(4.48) lim / fin{x,a)4>(x)dx= -{4>{a-)^^{a^)\. 

Since e in (4.10) and (4.11) can be arbitrarily small, (4.8) now follows from (4.48). 
This completes the proof of the theorem. □ 
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5. Airy and parabolic cylinder functions 

We now turn our attention to the integral representation (1.7). Here, the interval 
of concern is the whole real line. However, the argument for this result remains 
similar to that for Theorems 1 & 2, and we will keep it brief. As before, we let b 
and r] be positive numbers such that b > max{l, \a\}. The Airy function Ai(<; — x) 
is a solution of the equation 

(5.1) ^ + {x-t)y = 0, -oo<i<oo, 
and satisfies 

(5.2) Ai(oo) = Ai'(oo) = 0, 

(5.3) Ai{-t~x) = -^{t + xr'^/'^[sin({x,t) + ei{x,t)], 
and 

(5.4) Ai'{-t-x) = — ]={x + ty/^[cosC{x,t)+S2{x,t)], 



/TT 



where ei{x,t) = 0{t ^/^) and S2{x,t) = 0{t ^/^), as t ^ oo, uniformly for x > —b, 
and where 

(5.5) t:{x,t) = l{t + xr/' + '^. 

Define 

/oo 
Ai{t - x)Ai{t - a)dt. 
-R 

Prom (5.2), we have 

, Ai(-R-x)Ai'{-R-a)-AU-R-a)Ai'(-R-x) 

5.7 5R{x,a) = ^ '- ^ '- ^ ^ 

x — a 

see (3.6). 

Theorem 3. For any a € K and any piecewise continuously differentiable function 

(/>(x) on (—00,00), we have 

(t){x) / Ai{t - x)Ai{t - a)dtdx = -[(p{a~) + (/)(a+)], 
-00 J-R 2 

provided that the two integrals 



/ — I l-C 
\x\-i\(j){x)\dx and / 
-00 Ji 



X ^(j){x)\dx 



are convergent. 



Proof. Using the asymptotic formulas (5.3) and (5.4), one can show from (5.7) that 
there are positive constants Mi and M2 such that 



Jb 

and 



/ \SR{x,a)(j){x)\dx < Ml / x~^\4>{x)\dx 

Jb Jb 

/ — b j* — b 

\5R{x,a)(l){x)\dx<M2 / |a:|"s |^(a;)|da; 
-CSO J —00 
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for R > 2\a\ + 1. On account of the convergence of the two integrals in (5.9), for 
any e > there is a constant c > b such that 



(5.10) 



\dfi{x,a)(l){x)\dx < e 



for all R > 2\a\ + 1. 

In view of the asymptotic formulas (5.3) and (5.4), equation (5.7) also gives 



Tr{x — a)(5i{(a;, a) 



(5.11) 



R + a 
R + x 

R + x 
R + a 



sin ({x, R) cos ({a, R) 



cos({x, R) sin ({a, R) + e{x,a; R), 



where e{x,a;R) ^ 0{R-^/^), as R -> oo, uniformly for a;,a G (—6,6). In a manner 
similar to (4.29), by using the generalized Riemann-Lebesgue lemma we have 



(5.12) 



lim 



5ji{x,a)(j}{x)dx ~ lim 



sin{C(x,i?)-C(a,i?)} 



4>{x)dx, 



for any piecewise continuously differentiable function 4i{x) and any 77 > 0. Since 
(^(a;, R) — C,{a, R) = ^/R{x — a) + 0{1/^/R) for large R and bounded x and a, we also 
have sin{(^(a;, R) — C(a, R)} = sm{\/R{x — a)) + 0{\/\/R) as i? — !• 00 for bounded 
x and o. Thus, equation (5.12) gives 



(5.13) 



lim 



Se.{x, a)(f){x)dx — lim 



a+ri 



sin y/R{x — a) 



(j)(x)dx. 



By Jordan's theorem on the Dirichlet kernel [ , p. 473], the value of the last limit is 
[0(a-) + (j){a+)]/2. The final result (5.8) now follows from (5.10) and (5.13). □ 

The parabolic cylinder function W{a, x) is a solution of Weber's equation 



(5.14) 

with boundary conditions 

f2ka 



(5.15) 
and 
(5.16) 

where ka 
(5.17) 



dx'^ 



cosC(a, x) + O 



a y = 



X — ?> 00, 



2 /I 

- — sinC(a,a;) + 0( - 

kaX [ \X 

s/l + e2"« - e^'' and 

C{a,x) = - alnx + ^argpQ 

see [ I , p. 693]. The derivative of this function has the behavior 

1 



X — >■ 00, 



(5.18) 
and 
(5.19) 



W'{a,x) = 
W'{a,-x) = 



X 

2k n 



sinC(a, x) + O 
cos ({a, x) + O 
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This function is also related to the parabolic cylinder function U{a,x) via the con- 
nection formulas [ ' , p. 693] 

(5.20) W{a, x) ^ i2ka)hi'"' Re {e'(3"*^ + 5'^)i7(za, xe-i""')}, 

(5.21) W{a,-x) = (2/fca)5eJ™Im {e'(5<^^+5'')[/(ia, xe-J")}, 

where a; > and <j>2 = argr(i + ia). From the integral representation [12, p. 208] 

(5.22) U{a,z)^ / e— Rea>--, 
one can easily show that 

|t/(za,xe'^")|< x-i 
|r(^ +ia)| 

Since 

' ' (cosh7ra)V-i 

and 

v/2fc^-e-5™, W^-2e5''° 

V "-a 

as d — >■ +00, it follows that for large positive a, there is a constant Afi such that 
\W{a,x)\ < Miei''''x-^, x > 0, 

\W{a, -x)\ < Miei'^'^x'^, x > 0. 

As a — )■ —00, we have 

|r(i +m)| - %/2^e5'^" = %/2^e~5^l°l 



and 



^ 72, W |- - %/2. 



Hence, there exists a constant M[ > such that for large negative a, 

\W{a,x)\ < M[e^''\''\x-K x > 0, 

\W{a,-x)\<M[ei''^''^x-^, x > 0. 

By using (5.20), (5.21) and (5.22), it can also be shown that there are positive 
constants Af2 and M2 such that for large positive a, 

\W'{a,x)\< Ahei'^'x^, x > 0, 

\W'{a,-x)\ < Ahe^'^'x^, x > 0, 

and for large negative a, 

|l¥'(o,a;)| < Afae^'^l'^lcc^, x > 0, 

|W^'(a,-x)| < M^ei'^l^lxs, a; > 0. 

We define 

(5.23) SR{a,b)^[ W{a,x)Wib,x)dx. 
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From (5.14), one can derive 



W{a,x)W'{b,x)-W'{a,x)W{b,x) 
OR(a,b) = 



b — a 

Using the asymptotic formulas (5.15) and (5.18), we obtain 

{b - a)dR{a, b) - (^Kh + -^)sm{C(a, R) - ({b, R)} + O 

By an argument similar to that for Theorem 3, one can establish the following 
result. 

Theorem 4. For any a € M. and any piecewise continuously differentiable function 
4>{x) on (— oo,+oo), we have 

/+00 pR 
(f>{b) / W{a,x)W{b,x)dxdb = TT^/l+l^[(l){a-)+(l>{a+)], 
-oo J~R 

provided that two integrals 

|0(^)|ei-N^ and / \cj){x)\ei^^ — 

-oo FI Jl ^ 

are convergent. 

6. Series representations 

To prove the representations in (1.8) — (l-H), we only need to recall some ex- 
pansion theorems concerning orthogonal polynomials (functions). For instance, in 
the case of Legendre polynomials Pn{x), Theorem 1 in [9, p. 55] can be stated in the 
following form. 

Theorem 5. Let f{x) be a piecewise continuously differentiable function in (—1, 1), 
and put 

(6.1) Sr,{t,x)=f2{''+^^Pk{t)Pk(x). 

If the integral 

r-l 

/•2 



rix)dx 
1 



is finite, then 



1 

(6.2) hm / 5r.{t,x)fit)dt= -[fix+) + fix-)]. 

The statement in (6.2) is equivalent to that in (1.8); i.e., the finite sum in (6.1) 
defines a delta sequence. In a similar manner, one can restate Theorems 2 and 3 in 
[ , p. 71 and p. 88] as follows. 

Theorem 6. Let f{x) be a piecewise continuously differentiable function in (— oo, oo), 
and put 

n 

(6.3) <5„(t,a;)=e-*'^— -=i/fe(i)ilfe(a;). 

fc=0 '^''^ 
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// the integral 

/C30 
-C30 

is finite, then 

(6.4) lim / S„{t,x)fit)dt^-[f{x+) + fix-)]. 

Theorem 7. Let f{x) be a piecewise continuously differentiable function in (0, oo), 
and put 

(6.5) 5r.{t,x) = e-H'^Y. ^^^^[^^^ Lt\t)L'^\x). 
If the integral 



e 



-h'^f{t)dt, a>-l, 



/o 



is finite, then 

pea 

(6.6) lim / 5,,{t,x)f{t)dt=-[f{x+) + f{x-)]. 
To demonstrate (1.11), we recall the Laplace series expansion 

(6.7) /(^I'-Z-O^E^^/ f{02,h)Pk{cosj)smd2d62d^2, 



where cos 7 = cos^i cos 6*2 — sin 6*1 sin 6*2 cos((/)i — 02 ) and Pk{x) is a Legendre poly- 
nomial; see [ii, p. 147]. By the addition formula [ , p. 797] 

(6.8) P.(cos7) = ^— J2 YkiiOuMY:iie2,cb2), 

l = -k 

we can rewrite (6.7) in the form 

(6.9) f{ei,(bi) = f2 [ [ f (02,(^2) sm02Y.Yki{ei, ^i)Y:iie2, (b2)de2d^2. 

For / g C([0, tt] X [— TT, tt]), the series on the right converges pointwise to the function 
on the left; see [(i, p. 344]. This result can be expressed as follows. 

Theorem 8. Let f {61, (pi) he a continuous function on [0, tt] x [— tt, tt], and put 

n k 

(6.10) 5„{9i,e2)Sn{(j>i,(l)2) ■■= sm92Y, J2 0i)>7;(^2>2). 



fc=0 l = -k 

Then, we have 



(6.11) /(0i,</)i) = lim / / Sniei,e2)Sn{(puh)f{02,<l>2)de2dc^2. 

Equation (6.11) is equivalent to saying that (5„(6'i, 02)<5n(0i, ^2) is a delta se- 
quence of d{6i — 6*2)^(01 — ^2) for 9i,92 G [0,7r] and 01,02 G [— 7r,7r]. The use of 
the identity 5 (cos ^i, cos ^2) — sing ^(^1 ^ ^2) gives (1.11); see [s, p. 49]. 
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